Formation of Plasmoid Chains in Magnetic Reconnect ion 
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A detailed numerical study of magnetic reconnect ion in resistive MHD for very large, previously 
inaccessible, Lundquist numbers (10 4 < S < 10 8 ) is reported. Large-aspect-ratio Sweet-Parker 
current sheets are shown to be unstable to super- Alfvenically fast formation of plasmoid (magnetic- 
island) chains. The plasmoid number scales as S 3 ^ 8 and the instability growth rate in the linear 
stage as S 1 ^ 4 , in agreement with the theory by Loureiro et al. [Phys. Plasmas 14, 100703 (2007)]. 
In the nonlinear regime, plasmoids continue to grow faster than they are ejected and completely 
disrupt the reconnection layer. These results suggest that high-Lundquist-number reconnection is 
inherently time-dependent and hence call for a substantial revision of the standard Sweet-Parker 
quasi- stationary picture for S > 10 4 . 
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Introduction. Magnetic reconnection is a fundamen- 
tal plasma process of rapid rearrangement of mag- 
netic field topology, accompanied by a violent release of 
magnetically-stored energy and its conversion into heat 
and into non-thermal particle energy. It is of crucial 
importance for numerous physical phenomena such as 
solar flares and coronal mass ejections [lL |2|, magnetic 
storms in the Earth's magnetosphere 0,13, 111, and saw- 
tooth crashes in tokamaks @, 0] • Reconnection times in 
these environments are observed to be very short, usu- 
ally only 10 to 100 times longer that the global Alfven 
transit time, ta = L/va, where L is the characteristic 
system size and va is the Alfven speed. This is in direct 
contradiction with the classical Sweet-Parker (SP) @, 0] 
reconnection model, which employs the simplest possible 
non-ideal description of plasma — two-dimensional (2D) 
resistive magnetohydrodynamics (MHD) — and predicts 
a very long reconnection time scale t rec 

~ taS 1/i , where 

S = Lva/v ^> 1 is the Lundquist number, and n is the re- 
sistivity (or magnetic diffusivity) of the plasma. Both nu- 
merical simulations [l(J EH and laboratory experiments 
[T^ have confirmed the SP theory for collisional plas- 
mas where resistive MHD with smoothly varying (e.g., 
Spitzer) resistivity must apply. Because of this discrep- 
ancy between the MHD picture and observations, efforts 
to understand magnetic reconnection have moved beyond 
simple resistive MHD to increasingly sophisticated and 
realistic plasma-physics frameworks incorporating colli- 
sionless processes such as anomalous resistivity or two- 
fluid effects, where, indeed, fast reconnection rates have 
been found [HEi. For these reasons, simple resist ive- 
MHD reconnection has come to be viewed as well under- 
stood, uninteresting, and mostly irrelevant. 

However, most previous numerical studies of resistive 
MHD reconnection have been limited by resolution con- 



straints to relatively modest Lundquist numbers (S ~ 
10 4 ). The same is true for dedicated reconnection ex- 
periments. On the other hand, in most real applications 
of reconnection, the Lundquist numbers are much larger 
(e.g., S ~ 10 12 in the solar corona, S ~ 10 8 in large toka- 
maks). Thus, there is a large gap between the extreme 
parameter regime that we would like to understand and 
that accessible to the numerical and experimental studies 
performed so far. Asymptotically high values of S have 
never been probed by numerical simulations and, there- 
fore, one cannot really claim a complete understanding 
of magnetic reconnection even in the simplest framework 
of 2D MHD with a (quasi-) uniform resistivity. 

Of particular interest is the possibility that current 
sheets with large aspect ratios L/Ssp ~ where Ssp 

is the width of the SP current layer, should be unstable 
and break up into chains of secondary magnetic islands, 
or plasmoids — a phenomenon absent from the SP the- 
ory. A tearing instability of large-aspect-ratio current 
sheets was anticipated by Bulanov et al. [15] and Biskamp 
[lij. Current sheets were, indeed, found to be unstable in 
those numerical experiments where S ~ 10 4 was reached 
(e.g. 0, El, E3, m El HI). Current-sheet instabil- 



ity and plasmoid formation have also been observed in 
numerical reconnection studies using other physical de- 
scriptions, e.g., fully kinetic simulations I21L fijj [2^ , and 
there is tentative observational evidence [24|,[25| that they 
might play a key role in the dynamics of magnetic recon- 
nection in the Earth magnetosphere and in solar flares. 
In fusion devices, plasmoid formation is less well diag- 
nosed but recent results from the TEXTOR tokamak [26[ 
suggest that they might also be present. Theoretically, 
plasmoid formation has been proposed as a mechanism 
of fast reconnection 19|, 2tJ and non-thermal particle ac- 
celeration in reconnection events [23|. Thus, plasmoids 
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seem to be as ubiquitous as magnetic reconnection it- 
self. However, even though their appearance has been 
reported by many authors, neither the plasmoid forma- 
tion in the limit of asymptotically large S nor its effect 
on the reconnection process have been systematically in- 
vestigated on any quantitative level and remains poorly 
understood. 

As the first step towards this goal, Loureiro et al. 
[1^ developed a linear theory of the instability of large- 
aspect-ratio current sheets that, unlike in the calculation 
of Ref . [HI , emerges from a controlled asymptotic expan- 
sion in large S. Mathematically, the instability resembles 
a tearing instability with large A', leading to the forma- 
tion of an inner layer with the width Si nner ~ S~ 1 ^ 8 5sp- 
The instability is super- Alfvenically fast, with the max- 
imum growth rate scaling as jta ~ S 1 ^] the fastest- 
growing mode occurs on a scale that is small compared 
to the length of the current sheet, viz., the number of 
plasmoids formed along the sheet scales as S 3 ^ 8 . 

In this Letter, we report the next logical step towards 
the detailed assessment of the role of plasmoids in mag- 
netic reconnection: the first numerical evidence that not 
only do current sheets go unstable but that they do so in 
the extremely fast way predicted by Ref. [28[ and the in- 
stability quantitatively obeys the scalings derived there. 
To this end, we perform a set of 2D MHD simulations 
of an SP reconnection layer with uniform resistivity and 
asymptotically large Lundquist numbers 10 4 < S < 10 8 . 

Numerical Set Up. Probing such previously unattain- 
able values of the Lundquist number is made possible 
by a special numerical set up that effectively zooms in 
on the SP current sheet by choosing a simulations box 
whose size in the direction across the reconnection layer 
(x) is just somewhat larger than, but generally tied to, 
the SP thickness, L x > Jsp, while in the direction along 
the layer (y), the box covers a finite fraction of the global 
length L of the current sheet. The boundary conditions 
are used to mimic the asymptotic matching between the 
global and local solutions (in the spirit of 11]). Let us 
explain how this is done. 

We solve the standard set of compressible resistive 
MHD equations (the adiabatic index is 5/3; viscosity 
and thermal conductivity are ignored) in an elongated 
2D box, f 

Lx-> L x ] x [ Ly,Ly\. At the upstream bound- 
aries (x = dzijc), we prescribe the reconnecting com- 
ponent of the magnetic field, B y (x = ±L x ,y) = ±£>i n 
and the incoming velocity, v x (x = ±L x ,y) = =F^in- As 
the box is understood to model an SP current sheet, 
we set L x = Ssp = LS -1 ^ 2 = (Ltj/va) 1 ^ 2 , where va 
is the Alfven speed corresponding to B in and L is the 
(half-) length of the current sheet. We should then have 
^in = vaS -1 / 2 = (vaV/ 'L) 1 / 2 . We choose our code units 
so that va = 1 and L = 1. Then setting L x = r} 1 ! 2 and 
^in = V 1 ^ 2 enforces a fixed SP reconnection rate based on 
va — 1 and L = 1. Choosing L y = 1 would correspond 
to simulating the entire length of the current sheet, but 



it is clear that in this local set up only an inner fraction 
of the sheet can be computed accurately, so we choose 
L y = 0.24. At the downstream boundaries y = ±L y , free 
outflow boundary conditions are imposed. The method of 
characteristics is used to determine the remaining bound- 
ary conditions. Compressibility effects are minimized by 
ensuring that the Mach number M = va/c s -> where c s 
is the sound speed, is small (in our simulations, it never 
exceeds 0.1). The initial conditions are chosen so as to 
represent qualitatively an SP-like current layer (using the 
Harris sheet profiles). We do not choose an initial pertur- 
bation with a particular wave number; instead the system 
itself is allowed to pick the most unstable wave number. 

Time Evolution of the Instability. For S < 10 4 , the 
current density at x = settles down to a quasi-steady 
state, and no plasmoids are observed, consistent with SP 
theory. As the Lundquist number is increased, this pic- 
ture changes dramatically. The system does not settle 
into a steady state — instead, as predicted by the linear 
theory [28j], the layer becomes unstable and secondary 
islands (plasmoids) form, with reconnection occurring at 
multiple X-points. This is illustrated in Fig. [IJ where 
the time evolution of the current density for S = 10 7 
is shown. We see that a plasmoid chain develops along 
the sheet and also that the plasmoids closer to the cen- 
ter of the sheet grow faster than those farther away from 
it. In Fig. 02 we show the time evolution of the width 
of the plasmoid closest to the center of the sheet. At 
early times, it grows exponentially. The growth rates for 
different values of S are plotted in Fig. [H The scaling 
jta ~ S 1 / 4 predicted by the linear theory [28] is obeyed 
extremely well. The fact that the growth rate for the off- 
center plasmoids is a decreasing function of the distance 
along the sheet was not captured in the simple equilib- 
rium model used in Ref. [28] but it does emerge in the 
calculation for a general SP equilibrium [2^ . 

The linear stage ends when the plasmoid width ex- 
ceeds the width of the inner layer, S 

inner ^ S~ 1/8 S S P (see 
Fig. [3]). The subsequent nonlinear growth is slower, but 
still sufficiently rapid to permit the plasmoids to reach 
the width of the current sheet before being advected out 
by the mean outflow along the SP sheet, v y ~ (y/L)vA- 
That this outflow causes the plasmoids to drift outwards 
is illustrated by Fig. where we plot the time traces of 
the plasmoid O points along the sheet. The plasmoids 
that are further away from the center of the sheet move 
faster due to faster outflow. Note that new O-point traces 
appear in between already existing ones — we interpret 
this as evidence of secondary-plasmoid generation, aris- 
ing from the collapse of the X points between the primary 
plasmoids into secondary current sheets, which then go 
unstable (cf. [18]). The breaking of the secondary sheets 
is also visible in the last two frames of Fig. [U In some 
cases, we have also observed coalescence of nearby plas- 
moids (not shown). 

Once the plasmoids grow enough to touch the wall, 



FIG. 1: Contour plots of the current density showing the time evolution of an SP current sheet for S — 10 7 . The times shown 
are, from top to bottom, t = 0.63ta, t = 0.96ta, t = 1.09ta and t = 1.27ta- The domain shown is —5sp < x < 5sp (inflow 
direction, vertical), and — 0.12L < y < 0.12L (outflow direction, horizontal), where 5s~p — 3 x 10 4 is the SP layer width and 
L = 1 is the (half-) length of the current sheet (see text; only the central half of the simulation box is shown). 




FIG. 2: Contour plots of the current density for (from top to bottom) S = 10 4 (t = 4.5ta), S = 10 5 (t = 2.9r A ) and S = 10 6 
(t = 2.6ta); the case S — 10 7 is shown in Fig. [1] 



the simulation becomes invalid and is stopped. The fact 
that they do touch the wall means that in a more general 
global set up, they would grow to exceed the SP width 
#sp, so the current sheet is broken up. 

Spatial Structure of the Plasmoid Chain. The devel- 
opment of multiple magnetic islands for different values 
of S is illustrated in Fig. [2j an increasing number of 
plasmoids is observed as the Lundquist number increases. 
Fig.[6]shows the number of plasmoids vs. S. Again, there 
is excellent agreement with the scaling ~ S 3 ^ 8 of the most 
unstable wave number predicted by the linear theory [28[ . 
Note that to avoid possible boundary-condition effects, 
we only count the number of plasmoids present in the cen- 
tral half of the simulation domain, — 0.12L < y < 0.12L. 

Note that the decrease in the linear growth rate with 



y means that in practice the simulation must be run into 
a nonlinear state to observe the formation of the entire 
plasmoid chain. However, by the time the outer plas- 
moids become detectable (but still linear), the more cen- 
tral plasmoids can already be well into their nonlinear 
evolution, so some secondary plasmoid generation (and, 
in some cases, coalescence) might have already taken 
place. This means that there is a degree of imprecision in 
measuring the number of the primary plasmoids result- 
ing from the linear instability of the entire current sheet 
(the diagnostic in Fig. [6] counts the maximum number of 
plasmoids during the run), especially at the largest val- 
ues of S. However, the errors this introduces in Fig. [6] are 
not large and do not affect the validity of our claim that 
the theoretical scaling S 3 ^ 8 is obeyed. At even larger 
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FIG. 3: Time evolution of the half- width of the plasmoid 
closest to the center of the sheet for S — 10 7 . Half- widths of 
the inner layer, dinner [28j], and of the SP sheet, 5sp = L x , are 
shown for reference. 



\(j K? 10 G 10 v 

Lundquist number, S 

FIG. 4: The growth rate of the plasmoid closest to the center 
of the sheet vs. S. The theoretical slope S 1 ^ 4 28] is shown 
for reference. For S > 10 7 , we could not calculate the growth 
rates accurately because the linear regime was too short lived. 
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FIG. 5: Plot of the position of the O points vs. time for 
S = 10 7 . Dotted lines mark the times at which the current 
sheet is shown in Fig. [1] 



Lundquist numbers S > 10 8 , the dynamics along the 
sheet is likely to be even messier, so a statistical descrip- 
tion of plasmoid chains may be required (cf. 0, 

Discussion. The key questions that remain to be an- 
swered are how large the plasmoids grow after they ex- 
ceed the width of the SP layer, and what is their im- 
pact on the overall reconnection rate. Addressing these 
questions requires global (or, at least, intermediate-scale) 
simulations capturing regions both interior and exterior 
to the current sheet. We could not computationally af- 
ford such simulations and simultaneously resolve the very 
large values of S required to diagnose the current-sheet 
instability in its asymptotic form. In our view, it was im- 
portant, before undertaking a global reconnection study, 
to understand the nature of the instability. The conclu- 
sion of this Letter is that the instability exists, is super- 
Alfvenically fast (cf. [13, 0]) and, in the limit of large S, 



quantitatively follows the linear theory of Ref. 
We stress that this is the first numerical study that has 
been able to make this statement and thus, in a sense, 
demystify the phenomenon of multiple-plasmoid genera- 
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FIG. 6: The maximum number of plasmoids in the central 
part of the sheet, — 0.12L < y < 0.12L, vs. S. The theoretical 
slope from linear theory, S 3 ^ 8 ~[28l], is shown for reference. 



tion. This conclusion puts further investigations of the 
plasmoid effect on magnetic reconnection on a firm the- 
oretical footing. 

The body of numerical evidence for plasmoids is now 
so overwhelming that there should remain little doubt of 
their importance in reconnection processes. It is clear 
that for sufficiently large systems, plasmoid-dominated 
reconnection layers are inevitable. Plasmoid formation 
and magnetic reconnection are thus inextricably linked 
and further progress in understanding reconnection in re- 
alistic systems necessarily requires a theory that takes the 
plasmoid dynamics into account. It also requires experi- 
mental evidence — however, present-day magnetic recon- 
nection experiments have not yet been able to observe 
plasmoid formation, most likely because of the moderate 
current sheet aspect ratios. The effect of plasmoids on 
reconnection is, in our view, one of the natural objects of 
emphasis for the next generation of dedicated magnetic 
reconnection experiments. 
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